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I. INTRODUCTION 

At the classical level, quantum chromodynamics (QCD) 
with massless fermions is chirally symmetric; that is 
to say, the QCD Lagrangian is invariant under sepa- 
rate global flavor rotations of the left- and right-handed 
quarks. For Nj massless quark flavors the chiral symme- 
try group is 



SUviNf) X SUA{Nf) X Uv{l) X C/^(l) 



(1) 



where the vector (axial vector) symmetry groups de- 
noted by a subscript V {A) rotate left- and right-handed 
fermions in the same (opposite) way. The chiral symme- 
try is a good starting point for the Nf = 3 light quark 
flavors u,d,s. In the quantum theory, the SUA{Nf) is 
spontaneously broken by the dynamical condensation of 
quarks, (qq) ^ 0, which results in the generation of a con- 
stituent quark mass and gives rise to iV? — 1 Goldstone 
bosons, which can be identifled with the octet of light 
pseudoscalar mesons. Further, the small but finite (cur- 
rent) quark masses break the axial SU a{J^ f) explicitly 
and induce a mass for the pseudoscalar mesons. Further- 
more, the SUv{Nf) is softly broken by the differences in 
the current quark masses, which lifts the mass degener- 
acy of the pseudoscalar mesons. The t/y(l) symmetry 
corresponds to baryon number conservation and remains 
unbroken in the QCD vacuum. The Ua{^) symmetry, 
however, has been an issue for quite some time. If the 
t/yi(l) were intact, the light hadrons would occur in de- 
generate parity doublets, which is not the case. Further- 
more, if Ua{^) were spontaneously broken, there should 
be a (nearly) massless (or at least light) pseudoscalar fla- 
vor singlet meson, the "would-be" Goldstone boson of 
spontaneous UaI^) symmetry breaking. The only can- 
didate is the 77', which, however, is by far too heavy to 
qualify for this particle. 

It was first shown by Adler pl. Bell and Jackiw [4I 
that the ^^4(1) is anomalously broken, i.e., broken by 
quantum effects. The anomalous breaking of the global 
1/^(1) symmetry manifests itself in the non-invariance of 
the fermionic integration measure ISl, |4| and results in 
the well-known axial anomaly, which for massless quarks 
reads 



Here, ji^{x) = ip{x)j^"f5'ip{x) is the axial current and 
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(3) 



327r2 ^'•' 
is the topological charge density in Minkowski space, with 






(4a) 
(4b) 



being the field strength tensor and its dual (/"''^ are the 
structure constants of the su(iVc) algebra and g is the 
coupling constant; we use color sub- and superscripts 
indiscriminately). Adler and Bardeen showed |5| that 
Eq. ([2|) does not recieve corrections from higher-order di- 
agrams. Although q{x) is a total divergence 



F;,F^ = 2a^x^ 



where 



A:dpA: + -gr'^A';,A'pA^, 



(5) 



(6) 



is the topological current, there are topologically non- 
trivial field configurations (in Euclidean space) such as 
instantons [a, 0| , magnetic monopoles @ or center vor- 
tices 0, [l3| I for which the topological charge 



Q = d*x_E q{x) 



(7) 



is non-zero (for smooth field configurations of finite Eu- 
clidean action, Q is integer- valued) . As a consequence 
of the existence of these field configurations, the axial 
charge 



Q5 



d=^^J5°(^) 



(8) 



is not conserved [ll|, [l2| ■ 

Using large- iVc arguments, Witten jlj| and Veneziano 
|14| showed that the axial anomaly provides a mass term 
for the pseudoscalar flavor singlet meson given by 



d^j^ix)^ 2 Nfq{x). 



(2) 



rri^ — 2m^ 
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where F^^ — 93 MeV is the pion decay constant and are pure gauge spatial fields 

X^ f d^XE{0\q{x)q{0)\0) (10) A^AIta^ 



Ud,U\ 



9 



(12) 



is the topological susceptibility, which is a purely glu- 
onic quantity (in Minkowski space Eq. PT!]) has an addi- 
tional factor —i). Since this quantity is defined as vac- 
uum expectation value it is clear that its evaluation re- 
quires nonperturbative methods. Indeed, in perturbation 
theory the topological susceptibility vanishes to all or- 
ders. This quantity has been calculated on the lattice (see 
Refs. [la, [la] for recent calculations) and the results are 
compatible with the prediction of the Witten-Veneziano 
formula Eq. (O 



X ci (180 MeV)"^ 



(11) 



using the experimental data for the meson masses and 
Ftt as input. 

Like the string tension [l7[ , the topological susceptibility 
seems to be dominated by center vortices [l3|- In fact, a 
center vortex model of the infrared sector of Yang-Mills 
theory \\m yields a value for x compatible with lattice 
results [20| . 

Topologically nontrivial Euclidean field configurations 
such as instantons and center vortices describe quan- 
tum tunneling between topologically different Yang-Mills 
vacua [21|, see Eq. (fT^ below. Like in quantum me- 
chanics, this quantum tunneling is fully accounted for 
by the solution of the Yang-Mills Schrodinger equa- 
tion. Recently, progress has been made in determining 
the vacuum wave functional by a variational solution of 
the Yaiig-Mills Schrddin ger equ ation in Coulomb gauge 
[13, [H [3, m, [H, [13, mM 13 . There is good evidence 
to believe that the obtained wave functional contains the 
essential infrared physics: the absence of gluons from the 
physical spectrum in the infrared [25|, a linearly rising 
potential for static color charges [27[ and a perimeter 
law for the 't Hooft loop [28]. In the present paper we 
use the Yang-Mills wave functional determined in refs. 
[23 . [27| to calculate the topological susceptibility given 
by Eq. (dni). 

The organization of the paper is as follows: In the next 
section we briefly review the ^-vacuum in the Hamilto- 
nian approach. In section 3 we derive the expressions 
for the topological susceptibility in the Hamiltonian ap- 
proach in Coulomb gauge and evaluate the relevant ma- 
trix elements. Our numerical results are presented in 
section 4. Some concluding remarks are given in section 
5. 



II. THE e- VACUUM IN THE CANONICAL 
QUANTIZATION APPROACH 

Consider Yang-Mills theory in the Weyl gauge Ag = 0. 
The classical (time-independent) vacuum configurations 



where ta are the hermitean generators of the gauge group. 
Imposing the usual boundary condition that the gauge 
function [/(x) € SU{Nc) approaches a unique value for 
|x| —> oo (independent of the direction x) compactifies 
M.^ to 5*3 and consequently the C/(x) can be classified 
according to the winding number n[U] G n3(S'3) [2ll |. 
The classical vacuum configurations Eq. p^ belonging 
to different winding numbers are separated by infinite 
potential barriers. In the quantum theory, tunnelling be- 
tween the different classical vacua occurs and in a semi- 
classical picture, the barrier penetration is described by 
instantons, (space-)time-dependent solutions of the clas- 
sical Euclidean Yang-Mills equation of motion, which in- 
terpolate between vacuum configurations differing in the 
winding number by ±1. We will not resort here to a semi- 
classical description but instead approximately solve the 
Schrodinger equation, which fully accounts for the quan- 
tum tunnelling. 

In Weyl gauge the physical coordinates are the spa- 
tial gauge fields Af(x) and the corresponding canon- 
ical conjugate momenta are given by the chromoelec- 
tric field Ef{x) = Fq^{x). In canonical quantization 
the electric field is promoted to the momentum opera- 
tor n°(a;) = —i5/5A'^{x) satisfying canonical commuta- 
tion relations, and the Yang-Mills Hamiltonian reads 



H = 



1 



d\ 



(n^(x)) V (Bf (x))^ 



where 



B? 
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d,Al + '-f 



abc Ab 






(13) 



(14) 



is the non-Abelian magnetic field. 

In Weyl gauge. Gauss' law is lost from the equations of 

motion and has to be imposed as a constraint on the wave 

functional 



7?fn?(x)*[A]^-5p:xt(x)*[A] 



Here, 



jjab ^ gabg^ ^ ^^ab ^ ^ab 



racb AC 



(15) 



(16) 



is the covariant derivative in the adjoint representation of 
the gauge group and yOg^t (x) denotes the "external" color 
charge density of the matter fields. The operator DIl on 
the left-hand side of Eq. ([T5|) is the generator of time- 
independent small gauge transformations (ri[U] = 0). In 
the absence of external charges, p^xti^) ~ O7 Gauss' 
law requires the wave functional to be invariant un- 
der small gauge transformations only, while under large 
gauge transformation it needs to be invariant only up to 
a phase ,31, 32] 



^e[^c^] ^e-'''"!^!*^^] 



(17) 



Here, is a free real parameter, which characterizes the 
vacuum wave functional and is called the vacuum angle. 
Since n[U] G Z, the wave functional of the ^-vacuum 
Eq. pT]) has the property 



^e+2AA] =*e[^] 



(18) 



which qualifies 9 as an angle variable. 6 is not known a 
priori and is not determined by the theory itself but must 
be fixed by experiment. One measurable effect of 6 would 
be a non- vanishing neutron electric dipole moment [33|. 
Current measurements [3^ restrict 6 to the extremely 
smah value l^] < lO-^". 

The transformation property Eq. ()17|) can be realized by 
the ansatz 



^g[A] 



'ieW[A] 



[A], 



(19) 



where (j)[A] is a gauge invariant wave functional, (t)[A'^] — 
(j)[A], and 



W[A] = 



9 



167r2 
167r2 



d^xK°{x) 



(20) 



-ijk 



d^a 



A^d.Al 



y_ i-abc A a Ah ac 



is the Chern-Simons action, which changes under gauge 
transformations by the winding number n\U] 



W[A^] ^W[A]+n[U] 



(21) 



The wave functional Eq. ([19)) does, however, not fulfill 
Eq. P^ . In the appendix we show how Eq. P^ has 
to be modified to promote 9 to an angle variable. We 
also show there that Eq. P^ is a correct wave func- 
tional when is restricted to [0, 27r). In the following, the 
cyclic property Eq. (|18p of the vacuum wave functional 
will be irrelevant, since we are anyway interested only in 
infinitesimally small 0, see Eq. (P^ below. Therefore, we 
can use the simpler wave functional Eq. p9| . 
The quantity of interest is the topological susceptibility, 
Eq. (jlOp . which in the Hamiltonian approach can be de- 
fined by [13 



Vx 



AHH), 



d6l2 



no quarks 



^-0 



(22) 



where V is the spatial volume and 



{H)e = {^e\H\^g) , {^e\^e) = l (23) 

is the expectation value of the Yang-Mills Hamiltonian 
in the 0- vacuum. To evaluate the ^-dependence of (-ff)e, 
the following identity will be useful 



Mf (x) 87r2 ""^ ^""^ 



(24) 



where Bf is the magnetic field Eq. p^ . Obviously, the 9- 
phase in Eq. (|19p can only contribute to the kinetic term 



of the Yang-Mills Hamiltonian Eq. ^. With Eq. ([24 
we find from Eq. fl! 



ni-^e[A] = e-^''^[^l ( H^ - ^ Bf ) cj,[A] 



87r2 



(25) 



Inserting this relation into Eq. p^ and using the Bianchi 
identity 



DfBtix) = 



(26) 



we find that the wave functional (j)[A\ satisfies the same 
Gauss' law as ^e[^] 



DfiiiM]^-9P:.tM] 



(27) 



To make contact with previous results obtained in the 
Hamiltonian approach to Yang-Mills theory, see Refs. 
[251 [27j , it is convenient to work with the 6'-independent 
wave functional (t>[A\ and absorb the ^-dependence into 
the Hamiltonian by defining 



{^e\H\^e) = {4>\He 
Using Eq. ^^ wc find 

2 



H„ 



1 



H? 



9 

87r2 



Bf 



(Bf) 



(28) 



(29) 



An alternative way to arrive at this Hamiltonian is to add 
the topological 6'-term directly to the original classical 
Lagrangian 






9 pa pfiv 
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(30) 



87r2 



E B 



The 0-term is a total derivative and does not contribute 
to the classical equation of motion. It does, however, 
change the canonical momentum from H" = FQ^ to 



n? = f;^, + 1 



87r2 



Bf 



(31) 



and after canonical quantization in Weyl gauge one finds 
again the Hamiltonian Hg (l^5[) . 

Instead of working (for p°^j(x) = 0) with gauge invariant 
wave functionals 0[A], it is more convenient to explicitly 
resolve Gauss' law Eq. ([?f|) by fixing the gauge. For this 
purpose Coulomb gauge diA^ = is particularly con- 
venient and will be used in the following. In Coulomb 
gauge, the gauge field is transversal A = A-^ but this is 
not true for the momentum operator. We split the mo- 
mentum operator into longitudinal and transversal parts 
H ^ H^-fHll, where H-^ = -i5/5A^. The latter satisfies 
the canonical commutation relation for transversal fields 

[A^(x),Hf (y)] =i5''^f,,(x)5(x-y) , (32) 

where ty (x) = 5ij — didj/d^ is the transversal projector. 



Gauss' law Eq. ([27]) can be solved for the longitudinal 
part n" in the standard fashion yielding 



nil 



^gd{-bd)-\p,^t+Pg) 



(33) 



where p'l — A^ Hf- is the color charge density of 

the gluons and {~Dd) is the Faddeev-Popov kernel in 
Coulomb gauge. With the aid of Eq. ([55)) . one derives 
from Eq. ((29l) the gauge fixed Hamiltonian of the 9- 
vacuum (by considering {(f>\Hg\(J3) and using integration 
by parts in the kinetic term). This yields^ 



He — Hq 



87r2 



Hi 



2 \ 



(34) 



where 



i/o = \J^'^ [j-'un^)JUfi^) + Bf{^) i?f (x)] + 

+ ^ I d^x d^y J^- V"(x) J F'^'lAJi^, y)p'{y) 

(35) 

is the usual Coulomb gauge fixed Hamiltonian [35] for 
9 = and the ^-dependent terms are given by 

Hi - -^/d'^: [Sf(x)n,"(x) + J-in,"(x)^i?f(x)] + 

+ 1 Jd^xdh {G^'[A]i^,y)d!B:i^)gp\y)+ 

+ J-'gp'^{^)JG^'[A]i^,y)dfB^iy)] 



and 



H^ = ]^jd^xBt{^)Bt{^) 



(36) 



(37) 



In the above expressions, J = Det(— 9^1).;) is the 
Faddeev-Popov determinant, G is the Green's function 
of the Faddeev-Popov operator 

- d^Dfi^) G'^[A]i^,y) = <5"^ ,5(x - y) , (38) 

F denotes the Coulomb operator 



F'"'[A](x,y) = \{~d.Dr\-d'){-d,D,)-' 



ah 



(39) 



and p is the sum of external and dynamical color charge 
densities 



ahc A b -rjc 



a a , a a , rabc yi b tt' 

P = Pcxt +Pq^ Pcxt + / A n. 



(40) 



For the present purpose, the evaluation of the topologi- 
cal susceptibility x (which is entirely defined in the gluon 



^ In the following, all field and momentum operators are transver- 
sal and we will omit the symbol ±. 



sector), the external charges are not needed and we will 
put Pp^t = in the following. The Faddeev-Popov deter- 
minant J^ represents the Jacobian of the transformation 
from the (flat) non gauge-fixed configuration space to the 
(curved) space of the Coulomb gauge fixed fields. In par- 
ticular, this Jacobian enters the integration measure of 
the scalar product of wave functions 



($i,$2)= [vAJ^liA] ^2[A] 



(41) 



The integration of transversal field configurations ex- 
tends over the first Gribov region fi [30|, allowing the 
surface terms to be discarded in integration by parts. 
Although the integration should be restricted to the fun- 
damental modular region A C Jl [37|, there is evidence 
that integration over fi yields the same expectation val- 
ues [Hi. 

Following Ref. [251 we introduce the "radial" wave func- 
tional 



$[A] = J^/^ $[A] 



(42) 



which removes the Faddeev-Popov determinant from the 
integration measure of the scalar product. Matrix ele- 
ments of observables 0[A, 11] can then be expressed as 



/ PAJ $J O $2 = VA<^10^2 , 



(43) 



where we have introduced the transformed operator 

d[A, n] == ji/2 0[A, U]J-^^^ = 0[A, U] . (44) 

An operator 0[A] depending only on the field variable 
A^(x) is obviously not changed by this transformation, 
while the momentum operator transforms as 



n?f 



n^(x) 



i SlnJ 
2 M-(x) 



(45) 



The transformed Hamilton operator Hg = J^I^HgJ ^/^ 
is obtained from Bq by replacing H by H. Bq was ex- 
plicitly given in Ref. [25|. Furthermore, since B^iyi) is 
a function of the field variable only, H^ does not change 
{H2 = H-i). The explicit expression for Bi is given in 
the next section. 



III. MATRIX ELEMENTS FOR THE 
TOPOLOGICAL SUSCEPTIBILITY 

We are interested here in the topological susceptibility, 
Eq. (|22p . Since this quantity is defined as second deriva- 
tive of (Hq) at = 0, it is sufficient to calculate (He) up 
to second order in 9. This requires to treat Bi in first 
order and Bi up to second order perturbation theory in 
9. Wc then find 



Vx 



9 



H2 10) 



E 

n#0 



l(0|gi|^) 

En 



(46) 



where {\n)} denotes the set of eigenstates of Hq. We use 
here a Dirac notation for the radial wave functionals (|^^ . 
{A\0) = ^o[A]. 

For the unperturbed Hamiltonian Hq, we use the varia- 
tional results obtained in Refs. [2a, 123] ■ For the (radial) 
vacuum wave functional the following Gaussian form was 
chosen 



MA] - im 



— A/'exp 



d'xd'yAU^)co,,{^,y)A'^{y) 



(47) 



where u>ij{:x.,y) = iij(x)a;(x, y) and aj(x, y) is an inte- 
gration kernel determined by minimization of the energy 
density. By means of Wick's theorem we can express ex- 
pectation values of powers of field operators by the gluon 
propagator 



(0|A^(x)4(y)|0) 



1 



:ab , 



,'(x,y) 



(48) 



Note that a;(x, y) depends only on |x— y| and by Eq. ((48|) 

its Fourier transform represents the single quasi-gluon 

energy. 

The vacuum wave functional |0) Eq. ((47|l is annihilated 

af{x) |0) = by the operator 



1 

71 



ll^A'i + iujr^^^^m 



(49) 



which is the annihilation operator of a quasi-gluon with 
energy uj. Here and in the following we suppress the ex- 
plicit spatial dependence of the involved quantities and 
include the spatial coordinates in the Lorentz indices, so 
that contracted Lorentz indices imply integration over 
the spatial coordinate. The corresponding creation oper- 
ator reads 



.•^t 



1 

71 



ij 3 



■ -1/2-rTa 



(50) 



and from Eq. ([32]) follows that these operators satisfy the 
usual Bose commutation relation 

[<(x),af(y)]=<5'^''t,,(x)<5(x-y) (51) 

(temporarily restoring the explicit spatial dependence for 
clarity) . 

By repeated application of the creation operator a^ 
Eq. (|50p on the vacuum Eq. (I47p a complete basis for 
the gluon Fock-space of quasi-gluons is generated 



:AA„n<:^(x.)io>. 



(52) 



k=l 



For a proper normalization of these quasi-particle states 
the additional normalization factor J\fn is required when 
two or more indices take the same value, which can oc- 
cur due to the bosonic character of these quasi-particle 
excitations. 



With a complete basis of the Yang-Mills Hilbert space at 
our dipsosal we can explicitly carry out the perturbative 
calculations in Eq. (|46|) . The expectation value of H2 can 
be straightforwardly evaluated, yielding 






^g^EAKzllv 



16 



d^k d^q 3-(k-q)2 



(27r)3 (27r)3 cj(k)tj(q) 



(53) 



The second term in Eq. (j46p contains a sum over an in- 
finite number of states with an arbitrary large number 
of quasi-gluons and we have to resort to some approx- 
imations. Motivated by the variational calculation for 
9 = [2a, l23| we restrict ourselves to terms involving up 
to two loops in the energy. Under this assumption, Hi 
(cf. Eq. ([SSI)) is given by 

Hi = - fd'xIin^)Bn^) + 



d'x d'y G^'[A]{^, y) 9f Bf (x) gAf{y) ^y) 



(54) 



and we can use the following factorization in the matrix 
elements of the second term of Eq. ([5^ 

{0\G[A]{dB){gAU)\n) 

^{0\G[A]\0){0\{dB)igAn)\n), (55) 

where (G'[A]) is the ghost propagator. With these ap- 
proximations, the relevant contributions to the second 
term in Eq. (|iS|) 



n^O 



En 



come from two quasi-gluon states, n = 2 

1_ ^ |(0|i3n|2) - (G) (0| dBgAIl |2) p 



M2 = - 



2! ^^ E2 

1,2 ^ 



, (57) 



and from the 3-quasi-gluon states, n = 3 

Bn|3)p 



M, 



-y 



1,2,3 



E. 



(58) 



The prefactors 1/2! and 1/3! take the normalization of 
the states (I52p into account and avoid multiple counting. 
E2 and E^ , respectively, denote the energies of the two- 
and three-quasi-gluon states and accordingly are given 
by sums of, respectively, two and three w's. 
The matrix elements in Eqs. ([57)1 . ([55)1 can be evaluated 
by means of Wick's theorem, yielding 



M2 



N^ 



V 



d^fc 



g^M$^2lv 



(27r)3 w(k) 

d^fc d^q k^ 



M3 



2 

NcjN^, - 1) 
48 



V 



G(k + q) 



(27r)3 (27r)3 w(k) ' ^' \uj(k) Lu{q) 

d^fc d^q d^p w(k) + cj(q) + w(p) " 
(2^(2^(2^ w(k)t^(q)cj(p) 



1 + (k • q)2 k • q 



2 k^ ^ 

2(k.q)(k.p)(p.q)l (27r)3,5(k + q + p) 



(59) 



(60) 



r 



Within our approximations the expression in the bracket 
in Eq. ^ is given by the sum of Eqs. dSS]), §9^ and ((60)) . 
The term involving (G)^ in Eq. (j57p does not appear in 
Eq. ([551) since it involves three loops. It is worth noting 
that the first term in Eq. ((55)) cancels the first term in 
Eq. ([55]) . This cancellation is due to the fact that the 
Gaussian wave functional is exact in an Abelian theory 
and the two terms are actually the leading order contri- 
bution to the topological susceptibility in perturbation 
theory in powers of g, where x vanishes identically. In 
Ref. [331 it has been explicitly shown that this cancella- 
tion occurs also in next-to-leading order. 
The integrals in Eqs. ^, ([591) and §Q are UV diver- 
gent. Since the topological susceptibility x vanishes to 
any (finite) order perturbation theory in g, in principle, 
all UV-divergences should cancel. However, the approach 
[24 , |25| we are using is nonperturbative and, due to the 
approximations involved, does not include all terms of 
a given power in g. As a consequence, there is a mis- 



match of UV-singularities and our expression for x is 
UV-divergent. To remove these spurious UV-singularities 
we subtract from all propagators the corresponding tree- 
level form, i.e., we make the following replacements 



G(k) 



^(k) 



^(k)-i/Vk 



G,(k) = G(k) - l/k^ 



(61a) 
(61b) 



This is in the spirit of the zero-momentum subtraction 
scheme. The replacement (|6T|) makes the integrals (|53|) . 
([59]) and ([60]) convergent. Using the symmetry of the in- 
tegrands we can finally cast the expression for the topo- 
logical susceptibility into the form 



X = 9 



2 \ 2 
£^j N,{Nl^\){h+h). 



(62) 



where 



'^ = 47 (2^(2^"^ (''^"^ (^^ 



,_1^2l-(k-q) 



(k + q)2 



h^- / (03 (03fc'^.~'(k) k-^(k)+a;7'(q)]G.(k + q) 



1 + (k • q)2 q • k 



(63a) 
(63b) 



r 



Eq. ()63ap follows from the sum of Eq. ([60]) and the second 
term in Eq. ([55]) while Eq. (|63bp follows from the second 
term in Eq. ([59l) . The angular integration in Eq. (|63a|) 
can be performed analytically. 



IV. RESULTS 

The expressions for the topological susceptibility, 
Eqs. (|62l) and ((63)) . depend via the gluon and ghost prop- 
agators, cj^^ and G, on the vacuum properties, i.e., on 
the vacuum wave functional. We will use here the ghost 
and gluon propagators determined in Ref. [27[ as input. 



Furthermore, to simplify our calculations the numerical 
solutions obtained in Ref. [23| for the gluon and ghost 
propagators were fitted by the following ansatze 



c.(k) = \/k^ + ^ 






(64a) 
(64b) 



{g is the coupling constant). The parametrization (j64a[) 
of the gluon energy was already assumed heuristically by 
Gribov 36] . A factor g was included in the parametriza- 
tion of the ghost propagator (j64bp in order to facilitate 
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FIG. 1: Comparison between the numerical results of Ref. [23] 
(crosses) and the parametrizations (|64|) (hnes) for the inverse 
gluon propagator (top panel) and the ghost form factor (bot- 
tom panel). 



the fitting to the numerical results of Ref. [27|, where a 
factor g was included in the definition of the ghost form 
factor d(k) ^ 






(65) 



The numerical solutions of Ref. [27| and their parame- 
trizations by Eq. ([M|) are shown in Fig. [T] The integrals 
Eq. ([63)1 entering the topological susceptibility Eq. ([62)1 
recieve their dominant contributions from the infrared 
momentum regime, where the parametrizations Eq. (|64|) 
give a perfect fit to the numerical solutions. The infrared 
part of the parametrizations Eq. (|64p 



wiR(k) = 



rfiR(k) = 



M 



(66) 



was fitted to the infrared behaviour of the gluon and 
ghost propagators found in Ref. |27| by solving the cor- 
responding Dyson-Schwinger equations (DSEs). In the 
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FIG. 2: Running coupling as calculated in [27 



Hamiltonian approach in Coulomb gauge the physical 
scale is the Coulomb string tension ac^ i-e., the coeffi- 
cient of the linear term in the non-Abelian Coulomb po- 
tential. In units of the Coulomb string tension ac the fit 
of Eq. (|66p to the numerical solutions of the DSEs yields 
m^ = 0.614 and M = 4.97, while from the (analytic) in- 
frared analysis of the DSEs [2^ [23, |4l| one extracts the 
values m^ = 2/7r ~ 0.637 and M = \/8^ ~ 5.01. 
With the definition Eq. ([S5)) of the ghost form factor in 
the variational calculation of Refs. [22, [H [H, [H, [H, [l^l 
the coupling constant g drops out from the DSEs and the 
value of the coupling constant never had to be specified. 
Contrary to this, the topological susceptibility Eq. (|S^ 
explicitly contains the coupling constant. In principle, af- 
ter a complete renormalization procedure, in all physical 
(renormalized) quantities the coupling constant should 
be replaced by the running one defined in a renormaliza- 
tion group invariant way. Such a complete renormaliza- 
tion program is not yet feasible, although some progress 
in this direction has been made^. Fortunately the run- 
ning couling constant calculated in Ref. [23] in the Hamil- 
tonian approach in Coulomb gauge (see Fig. ^ has a very 
weak momentum dependence in the infrared regime. It 
basically stays constant below the infrared scale y^. It 
is this low-momentum regime which gives the major con- 
tribution to the integrals Eq. (|63p . We will therefore 
use the plateau value of the running coupling constant. 
Using the definition of the nonperturbative running cou- 
pling given in Ref. [43|, its value at zero momentum is 

El 



,n^ 5^(0) 167r 



(67) 



Numerical (Gauss-Legendre) evaluation of the integrals 



^ The inverse ghost form factor gd ^ (k) has the meaning of the 
dielectric function of the Yang— Mills vacuum [4011 . 



^ The counter terms required for the renormalization have been 
identified [13, [H- 



8 




FIG. 3: Value of X"^ /V^ depending on the ratio Ocjo. The 
two horizontal lines limit the range of the lattice results. 



Eq. ^ yields 

/i = (0.077 VCT;)^ 1-2. = (0.021V^)^ (68) 
Then we find for the topological susceptibility for Nc — 2 



0.45. 



(69) 



It was shown in Ref. [44!] that CTc is an upper bound for the 
string tension a extracted from the Wilson loop. Lattice 
calculations performed in Ref. [4^ for SU{2) show that 
ac — 1.5 (J, while Ref. [4y| seems to suggest an even larger 
value for a^ In Fig. [3] we present our numerical result for 
the topological susceptibility as a function of ac/(J- For 
1 < (Tc/cr < 1.33, Y is inside the range predicted by the 
lattice data [13, S, Hi, ^] . Choosing CTc = 1.5 a we find 
with y/a = 440 MeV 



X = (240 MeV)^ 



(70) 



This value is somewhat larger than the lattice prediction 
X^ (200 - 230 MeV)l 



reasonable agreement with the lattice data. The results 
obtained in the present paper are quite encouraging for 
the calculation of further hadronic quantities within the 
present approach. This will require the inclusion of the 
quarks. 
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APPENDIX 

In the following we show how the wave functional Eq. ^ 
can be modified to fulfill Eq. p^ . so that 9 becomes a 
true angle variable. For this purpose consider the wave 
functional 

pin 

^g[A]= / d9'd[e,9']^e'[A], (A.l) 

Jo 



where 



'•"■i-^S 



-im{e-e') 



(A.2) 



is the periodic (5-function satisfying d[9 + 2TT, 9'] = d[9, 9']. 
Inserting Eq. P^ into Eq. (|A.1[) we obtain 

MA] - ^e"'«™/(m - W[A]) cl,[A] , (A.3) 



V. SUMMARY AND CONCLUSIONS 

We have calculated the topological susceptibility x 
within the Hamiltonian approach to Yang-Mills theory 
in Coulomb gauge using the gluon and ghost propaga- 
tors determined previously by a variational solution of 
the Yang-Mills Schrodinger equation as input. To our 
knowldedge this is the first ab initio continuum calcu- 
lation of the topological susceptibility, granted the ap- 
proximations adopted. The precise numerical value for 
X depends on the relation between the Coulomb string 
tension and the one extracted from the Wilson loop. It 
is therefore very desirable to perform a sophisticated cal- 
culation of the Wilson loop within the present approach. 
Adopting the relation ac = 1.5 cr as suggested by lat- 
tice calculations, the numerical value obtained for x is in 



where 



fix) 



1 

2^ 



(A.4) 



The functional *e[^] satisfies both Eq. dTT]) and Eq. HH]) 
and thus can be used as the wave functional of the 9- 
vacuum. According to Eq. (|A.3p . it has the form 



§4A]-^e-*''™vl,^[A], 



(A.5) 



where 



<i'„,[A]^f{m-W[A])(f>[A] 



1 

2^ 







is a localized wave functional centered at the classical 
vacuum UdU^ with winding number n[U] = m. It satis- 
fies the relation 



*™[A^] =*„_„[[/] [A] 



(A.7) 



Replacing ^0[A] by \l'm[A] corresponds to the so-called 
"tight binding" approximation in solid states physics. 
We are not using this approximation. From Eq. (|A.6p 
it is seen that our wave functional ^e[A] is the Fourier 
transform of ^m [A] . 

Using the Poisson relation, the periodic (5-function 
Eq. (|A.2p can be expressed as 



6[e,e'] = J2 s{0-e' + 2Tri) , 



(A., 



l=-c 



where S{x) denotes the ordinary 5-function. Inserting 
this representation into Eq. (|A.ip we find 



^e[A]^^e+2.io[A], 



where Iq is defined by the condition 



+ 27r/o e [0, 27r) 



(A.9) 



(A.IO) 



This shows that vE'9[A], Eq. (fT9|) . is a correct repre- 
sentation of the wave functional of the ^-vacuum for 

0e [0,27r). 
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